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The generation and destruction of the supercurrent in a 
superconductor (5') between two resistive normal (A'^) current 
leads connected to a current source is computed from the 
source equation for the supercurrent density. This equation 
relates the gradient of the pair potential's phase to electron 
and hole wavepackets that create and destroy Cooper pairs in 
the N/ S interfaces. Total Andreev reflection and supercur- 
rent transmission of electrons and holes are coupled together 
by the phase rigidity of the non-bosonic Cooper-pair conden- 
sate. The calculations are illustrated by snapshots from a 
computer film. 

PACS: 74.25. Fy, 74.80.Dm, 74.80.Fp 



I. INTRODUCTION 



Andreev scattering (AS) of electrons into holes and 
vice versa by spatial variations of the superconducting 
pair potential in competition and cooperation with 
conventional scattering, determines the electronic struc- 
ture and transport properties of inhomogeneous super- 
conductors. The Tomasch effect in tunnel junctions 
1^, Josephson currents excess currents, and sub- 

harmonic gap structures [^pQ in superconducting {S)- 
normal conducting (iV)-superconducting junctions, as 
well as the transfer of half of the Magnus force to the core 
electrons of a moving vortex line are due to AS. It is 
involved in the persistent currents around the Aharonov- 
Bohm flux in an N/S metal loop iQ, and there is AS 
in He"^, too A wealth of AS phenomena has been 

discussed recently in and . 

While the conversion of a normal current into a su- 
percurrent by electron hole scattering in the interface 
between an N and an S region of semi-infinite lengths has 
been described before [ p^ , p^p |], the reverse process has 
not. Thus, it is the purpose of this paper to analyze the 
normal-current < — > supercurrent conversion processes in 
a superconducting layer of finite length Lz between two 
normal current leads. These normal leads are connected 
to a reservoir ( "battery" ) which acts as the current source 
in the closed circuit. The extensions of the N and S re- 
gions in X- and y-directions are and Ly] assuming 
that they don't exceed the London penetration length 
one may neglect inhomogeneities of the current density 
in the superconductor. The metal-vacuum boundaries 
are treated as rigid walls. By varying Lx and Ly one 
can vary the dimensionality of the system. By showing 



in detail how in any transport experiment involving su- 
perconductors electron ^ hole scattering brings about 
normal current supercurrent conversion our analysis 
may also prove useful for the understanding of trans- 
port phenomena in quasi-two-dimensional (Q2D) super- 
conducting/semiconducting heterojunctions |20-22| and 
superconducting quantum dots in Q2D channels. 



II. CHARGE CONSERVATION 

AS and the associated formation and destruction of 
Cooper pairs and supercurrents can be calculated from 
the Time-dependent Bogoliubov-de Gennes Equations 
(TdBdGE) They describe the evolution 

of the spinor quasiparticle (q.p.) wavefunction with 
the electron component M„(r, t) and the hole component 
D„(r, t) under the influence of scalar and vector potentials 
V{y. t) and A(r, t) in the single-electron Hamiltonian 
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-V - eA(r,i) 



+ Vir,t)-^l 



via the matrix equation 



(1) 



Here, the matrix Hamiltonian T-C{r,t) has Ho{r,t) and 
— iJQ(r,t) in the diagonal, and the pair potential A{r,t) 
and its complex-conjugate in the off-diagonal. The chem- 
ical potential /i in Hq is that of the reservoir. We neglect 
all influences of entropy production associated with cur- 
rent flow on the chemical potential, because the number 
of degrees of freedom of the reservoir is assumed to be 
much larger than that of the normal leads and the super- 
conductor. Thus, fj, is constant in space and time P3[ . 
The index n characterizes the stationary q.p. states from 
which the solutions of eq. (|l|) evolve after time-dependent 
scalar and vector potentials have been switched on. 

AS is a many-body process. For its analysis it is con- 
venient to consider a non-equilibrium conflguration \Ti^) 
of the many-body system where one quasiparticle state 
{la), characterized by a tripel I of quantum numbers and 
spin cr, is definitely occupied and all other q.p. states 
(na) are occupied according to the equilibrium distri- 
bution function /„. All interactions that might affect 



the spin are neglected. Then, it has been shown 23 
with the help of the TdBdGE (^ that the expectation 
values {Ti^\p{r,t)\Ti„) and {Ti„\i{r,t)\Ti„) of the many- 
body charge- and current-density operators satisfy the 
relation 
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divj. 



The electron and hole wavefunctions ui and vi satisfy 
eq. and ]si is the supercurrent density induced by 
the momentum and charge transfer from the q.p. in (la) 
to the superconducting condensate. All mean-field q.p. 
states in jsi and in the selfconsistency equation for A(r, t) 
are in a Hilbert space specified to \Ticr). Therefore, their 
wavefunctions all acquire the same phase shift Si{r,t) 
caused by the q.p. in {I, a). This leads to a phase shift 
2Si{r,t) of the pair potential, and j^; becomes propor- 
tional to the gradient of Si (r, t) times the number N of 
electrons in the superconductor |18 2^. (Recently the 
necessity of a phase gradient for charge conservation in 
7V/5'/A'^ junctions has again been pointed out by Sanchez- 
Cafiizares and Sols |2j].) The requirement that charge is 
conserved in the many-body system results in the funda- 
mental source equation for the supercurrent density 

dWU = 4|lm[A*(r,t)wi(r,t)z;r(r,t)](l - /,). (2) 

This equation has a non- vanishing right-hand 

side if the ui and vi describe quasiparticles that decay ex- 
ponentially in the superconductor during total Andreev 
reflection, because their energies are less than the maxi- 
mum value of A. In this case the source equation yields 
a finite supercurrent j^;. The phase shift of the pair po- 
tential 25;, on the other hand, is essentially given by the 
(integral of the) r.h.s. of eq. (||), divided by 1 
p8[. Thus, despite of its importance for charge conser- 
vation, its numerically tiny value can be neglected in our 
calculation of solutions of eq. (j^). 



III. CURRENT FLOW AND REPRESENTATIVE 
WAVEPACKETS 

A shifted Fermi sphere (or its equivalent in Q2D and 
QID conductors) represents the current-carrying many- 
body configurations in the two parts of the normal cur- 
rent leads that are parallel to the 2:-axis and connected 
to the superconductor. (The direction of current flow in 
the parts bent towards the reservoir is irrelevant.) These 
leads, supposed to be much longer than the mean in- 
elastic scattering length, are conductors with resistance 
psf . In this non-equilibrium distribution of electrons 
above the Fermi surface in states with positive momen- 
tum in z-direction and unoccupied states with negative 
z-momentum below the Fermi surface the current-driving 
force from the battery is balanced by the frictional forces 
from the energy-relaxation processes. The quasiparti- 
cles in this resistive-state configuration are uncorrelated. 
Thus, the total current in the closed circuit is the sum of 
the currents from the individual quasiparticles. 



In the following we try to obtain the details of nor- 
mal current < — > supercurrent conversion by studying 
the motion of the electrons (+) and holes (-) that are 
part of the shifted Fermi sphere. Their z-momenta are 
hk^{E), with k^{E) = [klp±E2m/h^Y/'^. Here, k.p = 
[fc|, — {tIxIt/LxY ~ {''^y''^ / LyYY^"^ is the z-component of 
the Fermi wavenumber kp = [2m/i/fi,^]^/^, and {uxTt/Lx) 



and (nyTr/Ly), 



integers, are the wavenumbers of the 



standing waves between the rigid walls that limit the 
metals in x- and y-directions. The energy E of both 
electrons and holes is positive and measured relative to 
the surface of the unshifted Fermi sphere at the chemi- 
cal potential fi. For normal current densities below the 
critical current densities of conventional superconductors 
all E are less than the modulus A of the pair potential 
A(z) w A • e(z)e(Lz - z); here e(z) is the Heavyside 
function which is sufficient to model the spatial varia- 
tion of the pair potential in the context of current flow 
[||-0,^4|j2§j2^. Details of the change of A(z) close to the 
interfaces because of the proximity effect matter little in 
the integral of eq. (|^) that yields the supercurrent. 

The motion of wavepackets shows best the dynamics of 
quantum-mechanical processes. Thus, similar to the use 
of representative, well-controlled, preformed wavepackets 
in the calculation of differential cross sections for con- 
ventional scattering, we investigate the current dynam- 
ics in a closed N/S/N circuit by studying the motion 
of wavepackets that are representative for the electronic 
conflguration in any transport experiment involving a su- 
perconductor connected to a current source by normal 
conductors. Conventional scattering processes are dis- 
regarded because we are only interested in AS in the 
N/ S interfaces. In principle, impurity scattering could 
be treated with the help of the scattering matrix formu- 
lation 1^,^. This is especially convenient for devices 
that, unlike the ones considered by us, involve only a 
small number of incoming and outgoing channels. Inter- 
face barriers that weaken Andreev scattering have been 
considered previously |^,^,^; the competition between 
conventional and Andreev scattering has been discussed 
in terms of the diagonal and off-diagonal forces associ- 
ated with broken symmetries and a motion picture 
of the wavepacket dynamics in such a case can be viewed 
in the Internet under the address given in the caption of 
Fig. 1. If conventional scattering is present, the Gaussian 
spectral function D{E), see below, has to be multiplied 
by the probability amplitude of transmission. There are 
only few scattering impurities and no interface barriers 
present in a Q2D electron gas in a modulation-doped 
system consisting of an InAs channel between an AlSb 
substrate and a superconducting Niobium layer that in- 
duces a pair potential in the electron gas via the prox- 
imity effect [ pO| . For such an experimental set-up our 
calculations apply exactly (within the limits of the An- 
dreev approximation) . For the general case of any super- 
conductor between any two normal current leads they 
show the quantum mechanically and electrodynamically 
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essential dynamics of AS that rules the charge transport 
in addition to conventional scattering. 

We start our analysis with the initial condition that a 
normalized electron wavepacket, localized around zq < 
in the normal current lead to the left of the superconduc- 
tor at time t ^ 0, travels towards the superconductor. In 
the center of the wavepacket the energy is Ei < A. (By 
varying Ei and kzp one can obtain all the low-lying elec- 
tron excitations that are part of the shifted Fermi sphere 
in the left normal lead.) For convenience we choose a 
Gaussian spectral function 



D{E) 



O-z ^ ^-\k+(E)-k+(Ei)?ai/2^-i\k+(E)-k+{Ei)Uo. 



the position- uncertainty parameter <IC |zo| is chosen so 
large that the related energy spread of the wavepacket, 
(5_E = h^k^{Ei)/maz, is less than A — Ei. Solutions of 
eq. (0) are calculated neglecting V and A and approx- 
imating A(r,t) by the real A(z), thereby disregarding 
repercussions of the quasiparticle-induced supercurrent 
on the q.p. jl^ and on itself. (As discussed above, the 
phase shift of the pair potential due to one Andreev re- 
flection process is negligibly small.) These solutions are 
multiplied by the spectral function D{E), integrated over 
all energies, and matched at the left N/S interface, i.e. 
z = 0, in the usual Andreev approximation where terms 
of the order of A/fi are neglected outside the exponen- 
tials. Energy-dependent functions are Taylor expanded 
around Ei up to first order in [E — Ei). This affects the 
amplitude of the Andreev-reflection probability 



— % arccos 



_E/A 



l{Ei)e 



{t/h){E-Ei)n 



and k^{E) = [k^p ± E2m/h'^]^/^ « kf ± {E - Ei)/hvf] 

Ti ~ h[A? — Ef]^'^/'^ is the time for one electron 
hole-scattering event and the associated formation of a 
Cooper pair, see eqs. (||,||), and 

kf = k^{Ei), = hkf/m. 

The resulting electron and hole wavepackets UNL{Y,t) 
and VNLi''c,t) in the left normal current lead, z < 0, 
and the exponentially decaying solutions us l (r, t) and 
^^^(rjt) in z > 0, that contribute to the source equation 
(0) essentially in the left half of the superconductor, turn 
out to be 

UNL = u;,e''='^"e-["''-"+"M'/2a^, (3) 
= u;i7(i;0e^'=r^e-[^«+("''/"r)^+",+ (*-^<)]V2a^.^ (4) 



usL = wie 
vsL = wa{Ei)e 



(5) 
(6) 



where 



Wl 
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iEit/h 



and Ki = (A^ - Ef)^^^ /Hvzf = I/WzF, with VzF = 
hkzF/'Ti- (For the sake of simplicity the complex 



wavenumbers in usl and vsl have not been Taylor ex- 
panded in {E — El) but rather taken at Ei right away.) 

Identifying the wavefunctions usl and from eqs. 
(5) and (H) with the ui and vi of the source equation 
(2) and integrating that equation from z = to z yields 
the density of the supercurrent in z-direction, induced in 
the left half of the superconductor by AS of the electron 
wavepacket unl into the hole wavepacket vnl, as 

jsLL ^ ez{2eVzF)\wi\^[l ~ e-^'^'^] 

y.g-{lzo+v^(t-r,/2)f + (TiV + /2f}/ali^^ _ j^y^ 



Here we have assumed that Lz ^ I/ki- In the opposite 
case one would have to add a second source term on the 
r.h.s. of eq. (^. This term would contain the contribu- 
tion from the wavepacket solution usR{'r,t),vsR{'r,t) of 
eq. (|l|) for < z < Lz that matches to the current- 
carrying q.p. wavepackets in the right normal current 
lead at z = and decays exponentially with increas- 
ing distance from Lz- However, in our case of large Lz 
the usi?Xr, t) and vsr(v, t) give only rise to the supercur- 
rent isi,R in the right half of the superconductor. The 
complex amplitudes of the solutions of the TdBdGE in 
the energy integrals that form the wavepackets zisij(r,t) 
and vsr(j, t) are uniquely determined by the requirement 
that the supercurrent densities jsi,L and jsi,Ri computed 
from usi^v*gj^ and usRVgj^, join smoothly at all times 
somewhere within the superconductor. Since the phases 
of the Fermi-liquid quasiparticles in the left and right nor- 
mal current leads are at random, only the current densi- 
ties, not the wavefunctions, must join smoothly. (If, on 
the other hand, the superconductor were only a thin slab 
with Lz <ti there would be a finite probability that 
the q.p. do not suffer AS and induce a supercurrent but 
rather carry their phases through the pair-potential wall 
in a tunneling process.) The matching point turns out to 
be Lz/2, and the evanescent wavepackets result to be 

usR = wie''''^'e-^'^^'-'^e-^''>+''^^'-'''^^^^^''% (8) 
VSR = wa{Ei)-^e''''^'e'^'^^^-''>e-^'"+''^'^^/^'''' . (9) 

More details are given in the Appendix. 

The supercurrent density js/.^, obtained by integrat- 
ing eq. (|^) from Lz to z, with usRVgj^ in the place of 
uiv'i , has the same form as jsi,L of eq. ^ except that 
exp[— 2k/z] is being replaced by exp[— 2k;(L^ — z)]. Fi- 
nally, the wavepacket solutions UNRiTc, t),VNRi'r, t) of eq. 
(0) in the right normal current lead, z > Lz, that match 
to the U5fl(r, t), vsr(jc, t) at the right interface in z = L^, 
become 



xe 



-lza+L^-z+v+ (t-Ti)f /2al 



-1 pik^pL^ iki (z-L^) 



VNR^wa{Ei) e"^""'e 
xe 



(10) 



(11) 
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Note that these wavepackets are the result of only one 
initial condition, namely "electron wavepacket incident 
from the left" . This initial condition, the matching of the 
wavepackets at the interfaces, the requirement of charge 
conservation as expressed by eq. (^, and the smooth join- 
ing of the supercurrent densities '^si,L and ]si,R determine 
unequivocally the wavepackets in the right normal cur- 
rent lead (apart from an irrelevant, constant phase factor 
that has deliberately been set equal to unity). 



IV. ANDREEV SCATTERING AND COOPER 
PAIR FORMATION 

Comparison of the unl, vnl from eqs. with the 

unr,vnr from eqs. ( po| , pd] ) shows that the center of the 
electron wavepacket unl, propagating to the right with 
velocity in the left normal current lead, and the center 
of the hole wavepacket vnr, propagating to the left with 
velocity in the right normal current lead, hit the left 
and right interfaces at z = and z — Lz (with different 
phases) at the same time — —zq/v'^, while the hole 
wavepacket vnl, propagating to the left in z < with 
v'l^ , and the electron wavepacket unr, propagating to the 
right in z > Lz, are retarded by the time t; with respect 
to the incident wavepackets. The holes moving to the left 
transport positive momentum and negative charge to the 
right, just as the electrons of opposite group velocity do. 
The supercurrent density j^; spreads "instantaneously" 
throughout the superconductor (the velocity of light not 
exceeding, of course) and couples directly electron hole 
scattering in the left to hole electron scattering in the 
right interface. 

This "instantaneous" coupling may seem to be surpris- 
ing, but only at a first look. Giving it a second thought 
one sees that our result confirms by explicit wavepacket 
calculation, and for the first time, as far as we know, what 
one has concluded before intuitively and from steady 
state calculations: hole electron scattering destroys 
and electron — > hole scattering creates supercurrents, 
whenever currents flow through S/N and N/ S interfaces, 
e.g. in vortex lines Both scattering processes must 
occur simultaneously even at far-apart interfaces because 
of charge conservation in closed circuits and as a direct 
consequence of the very essence of BCS superconductiv- 
ity: Cooper pairs are not bosons, despite of what one can 
read frequently, because within the volume of one Cooper 
pair of a conventional superconductor there are about one 
million of mass centers of other Cooper pairs, their cre- 
ation and destruction operators do not satisfy bosonic 
commutation relations, and their condensate wavefunc- 
tion is antisymmetric Therefore, the Cooper pair 

formed by total Andreev reflection in one interface can- 
not exist outside but only within the condensate to which 
the momentum 2hkzF and the charge of the two electrons 
of opposite spin, that have entered the superconductor, 
must be transferred. Because of the "phase rigidity of the 
electron pair fluid" |34|, so typical for off-diagonal long 



range order, this charge and momentum transfer, and 
the related phase shift, manifest themselves immediately 
in a current flow out of the other interface. Thus, no 
charges can accumulate in the superconductor according 
to its capacitance (as charges from quasiparticles with 
energies > A may do), and the reservoir has to emit 
holes into and receive electrons from the right normal 
lead at the same rate at which it emits electrons into and 
receives holes from the left current lead. This is illus- 
trated by the snapshots from a computer film in Figure 1. 
The Andreev-reflected wavepackets may be considered as 
supercurrent-transmitted wavepackets as well. Together 
with the incident wavepackets they represent any of the 
low lying single-particle excitations with energies Ei < IS. 
in the current-carrying, resistive parts of the circuit. 

Considering the supercurrent contribution from our 
representative quasiparticle-wavepacket configuration we 
note that its spatial maximum is at 2 = L^/l^ and its 
maximum in time occurs at t = to -f t;/2. From the 
equations for j^;, uslm. and vsl,r one sees that the quasi- 
particle current density 



jqpi = —Re 
m 



u, —Vui — vi — Vv, 
i i 



(1-/0 



changes into the supercurrent density jsi and vice versa 
within a distance 1/ki = Vzfti from the interfaces. 

In conclusion, the combinations of single-particle ex- 
citations (wavepackets) and collective modes (supercur- 
rents), connected by AS as shown in Fig. 1, are the con- 
sequence of the phase rigidity of the Cooper-pair con- 
densate and the adjustment of the current configuration 
to charge conservation in the N/S and S/N interfaces. 
They are the current-carrying elementary excitations in 
closed N/S/N circuits. 



V. OUTLOOK 

The supercurrent, carried by the condensate in the S 
layer, involves only states with |i?„| > A. It contin- 
ues the current in the normal current leads where all 
electrons and holes have energies Ei < A. If the total 
current density exceeds its critical value, i.e. if the cen- 
ter of the Fermi sphere in the normal current leads is 
shifted by more than hqcs = Am/hkp, depairing sets in, 
and when superconductivity has broken down, the un- 
correlated normal-state configuration reigns everywhere 
in the circuit. If, on the other hand, the single S layer is 
replaced by a mesoscopic SNS junction, the many-body 
configuration in the central N layer is a phase-coherent 
one and thus different from the uncorrelated configura- 
tions in the normal current leads. In an N/ SNS/N cir- 
cuit the SNS junction acts as a gapless superconduc- 
tor Q . It can carry a dissipation-free Josephson current 
through the central N layer via phasecoherent q.p. states 
with < A and \E\ > A This current converts as 
a whole into the total supercurrent of the S layers, and 
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vice versa, whereas, according to eq. (g), each uncor- 
related q.p. from the normal current leads individually 
induces its proper fraction of the total supercurrent. If 
the total current density exceeds the critical Josephson- 
current density at a Fermi-sphere shift of hqcSNS ~ h/d, 
where d is the length of the central N layer |5| , a voltage 
drop appears across the SNS junction. There are differ- 
ent models ||l^^,||,§ for SNS junctions with vohage 
drops. They differ with respect to the implicit assump- 
tions about the rate and energy range of q.p. creation 
in the central N layer by supercurrent — > quasiparticle- 
current conversion. The question of how this rate and 
range depend upon the weakening of phase coherence in 
the SNS junction by energy exchange between quasipar- 
ticles and electric fields, phonons, and thermal fluctua- 
tions like Nyquist- Johnson noise [^5| is presently investi- 
gated. 

APPENDIX: MATCHING OF SUPERCUR- 
RENTS 

The supercurrent density in the right part of the super- 
conductor, jsi,R, results from the evanescent wavepack- 
ets USR-, VSR- These are built up from the stationary 
solutions of eq. ([^) that decay exponentially within the 
superconductor with increasing distance from the right 
interface at z = L^. The free amplitudes n{E) of each of 
these solutions are determined by demanding that jsi,R. 
joins smoothly with the supercurrent density jsi.L from 
the left part of the superconductor at all times some- 
where within the superconductor. Once the amplitudes 
n{E) are known, the electron- and hole-wavepackets in 
the right normal current lead are unequivocally deter- 
mined by eq. and the matching conditions for the 
wavepackets in z = L^. 

The remainder of this appendix just shows how the 
amplitudes il,{E) are calculated. 

In the right part of the superconductor, with the Gaus- 
sian spectral function D{E) and those solutions of eq. (|^) 
that increase as z < Lz approaches Lz, we obtain 

usn = sin ( '^x\ sin ( '^y\ ^= ^~^E,t,n^^k.,. ^^.,iL.-.) 



Lx ) \ Ly J 5eV^^ 



Jo 

VSR = Sin (^x\ sin {'-^y] -1= e^^^'*/ V'=-^e--(^-^)7(i?0-' 
\ Lx J \ Ly / dEV2n 

-[{aJhv+Y{E~E, f/2+^{zo/hv++{t+Tl)/h){E-El)\ ^^^^ ^^^^ 



X / dEe 







(The term Se — h'^k^{Ei)/maz in the denominator, 
which results from replacing the wavcpackct integration 
over fc+ by one over E, drops out after the evaluation of 
the integrals.) 

Let z' be the point where the supercurrent densities 
jsi,L and jsi,R join smoothly. Then the integrals of the 
source equation (0) must satisfy 
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isl,L 



dz4-lm{A*usLV*sL){l- fi) 



dz 4-Im {A*usRV*sj^) (1 - //) = 3si,R 



Finally, comparison of the prefactors in eqs. ( |A5| ) and 
(0) yields 



(A3) 



(All) 



For convenience we write the complex n{E) as the ^hus, the complex amplitudes of the solutions in the 
product of two factors, one of them being energy- 
dependent: 



n{E) = ujiUj2{E). 



(A4) 



right part of the superconductor are given by 
2 



n{E) 



J{E~E,)T,/h 



Now we insert the usr and vsr of eqs. (Al) and ( |A2| ) 
into the supercurrent density js;^fl|z' of eq. (A3) and, sim- 
ilarly, jsi,L\z' is expressed by the solutions usl and vsl 
in the same form, i.e. without evaluating the energy in- 
tegrals: 



isLL 



-Asm"' 

h 



dEe 



-^{aJhv+y{E~Eif/2+i{zo/hv++t/h){E-Ei)\ 



dEe 



-[{a,/hv+Y(E-E,)''/2+i{zo/hv + + {t-Ti)/h)(E^Ei)j 



dze 



-2kiz 



(1-//), 



(A5) 







3sl,R 



-ez4-Asin^ 

h 



E. 



^1 



xlm<^ [7(^0- 



dEe 



-[{aJhv + Y(E-E,f/2+i{zo/hv ++t/h){E~Ei)\ ^^^^^ 



dEe 



-[{a:,/hv+YiE-Ei)'^/2+i(zo/hv + + {t+Ti)/h)iE-E,)j ^^^^^ 



dze-2-<(i--)(l_/;). 

z' 

Note that 

lm[[^{EO-'Y} = -lm{^{Eir}. (A7) 
We demand that the integrals over z are equal at z': 



(A6) 



dze--"^"- 



2ki 



1-e 



-2niz 



! e 



-2kiL^ 



2ki 



^2k[ z 



This equation is satisfied by 



1 



z' = —L,. 
2 



(A8) 



(A9) 



The energy integrals in eqs. ( |A5| ) and (A6) are real. By 
comparing them, i. e. the first one in eq. (A^) with the 
second one in eq. ([Aq) , or the second one in (A5) with 
the first one in eq. (|A6[), we find 



(A12) 



0J2{E) 



_ i{E-Ei)Ti/h 



(AlO) 
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FIG. 1. Propagation and Andreev scattering of the proba- 
bility densities \uf and |up of a representative spin-up elec- 
tron- and hole-wavepacket configuration (solid lines), and 
the induced supercurrrent density js (dashed line), in a cur- 
rent-carrying closed N/S/N circuit. For the sake of clear 
graphical representation the initial conditions for the elec- 
tron wavepacket, incident from the left, have been chosen 
as: energy in the center Ei = A/2 = 0.15 meV, spatial 
spread — 5 /xm, and k^p = 0.9fcF = 2.06 nm"'^. Via 
charge conservation by supercurrent induction these condi- 
tions determine the parameters and the motion of the re- 
sulting hole wavepacket incident from the right; the retar- 
dation time for AS, n = h/[A'^ — Ef]^''^, is about two pi- 
coseconds (ps). More computer films on electron — > hole and 
electron — > electron scattering in one N/S interface, also for 

energies above the gap, can be viewed in the Internet under 

http:/ /theorie.physik. uni-wuerzburg. de/TPl /kuemmel/filrns/filmse.htm, 
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